Spinodal dewetting of thin films with large 
oo ' 
O : interfacial slip: implications from the dispersion 

relation 

K ~ 5 ■ Markus Rauscher, a6 *Ralf Blossey, c Andreas Munch, d 

^ and Barbara Wagner 6 

i—i' a Max-Planck-Institut fur Metallforschung, 

<+h Heisenbergstr. 3, 70569 Stuttgart, Germany, 

b ITAP, Universitat Stuttgart, 
& . Pfaffenwaldring 57, 70569 Stuttgart, Germany, 

£h c Interdisciplinary Research Institute, c/o IEMN Avenue 

'"O ■ Poincare BP 60069, F-59652 Villeneuve d'Ascq, France 

Q | ^School of Mathematical Sciences, 

,^/ University of Nottingham, NG7 2RD, UK 

^ H \ e Weierstrass Institute for Applied Analysis and Stochastics, 

>- ■ Mohrenstrasse 39, 10117 Berlin, Germany 

^t" ! 

>5 ■ Version of July 16, 2008 

in ■ 
(N; 

o; 

00 . 

^^ ■ We compare the dispersion relations for spinodally dewetting thin liq- 

uid films for increasing magnitude of interfacial slip length in the lubri- 
. ,— ( , cation limit. While the shape of the dispersion relation, in particular the 

j^ ■ position of the maximum, are equal for no-slip up to moderate slip lengths, 

?H ' the position of the maximum shifts to much larger wavelengths for large 

slip lengths. Here, we discuss the implications of this fact for recently de- 
veloped methods to assess the disjoining pressure in spinodally unstable 
thin films by measuring the shape of the roughness power spectrum. For 
PS films on OTS covered Si wafers (with slip length b (a 1 /jxa) we predict 
a 20% shift of the position of the maximum of the power spectrum which 
should be detectable in experiments. 
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1 Introduction 

Wetting and dewetting phenomena are not only part of our everyday life 
but they are particularly relevant to technological applications (e.g., in 
coating processes) and in biological systems. The dynamics of films of 
a thickness smaller than 10 or 20 nanometers is not only governed by 
hydrodynamics, but the finite range of intermolecular forces, which are 
responsible for the richness of wetting phenomena [UH], becomes relevant 
[3J. This is true in particular for spinodally unstable films which have 
been analysed quantitatively in the framework of the thin-film equation 
[4]. In addition, the quantitative analysis of he roughness power spec- 
trum has been used in order to measure the disjoining pressure (DJP) (or 
the effective interface potential) between the liquid-solid and the liquid- 
vapor interface which is a result of the interplay beween the interactions 
among the fluid molecules and the interactions between the fluid and the 
substrate [5] [6] [7J. 

In the spinodally dewetting systems studied in Refs. [6] [7] [4], i- e -i 
polystyrene (PS) on silicon (Si) wafers covered with a native oxide layer, 
hydrodynamic slip between the fluid and the solid substrate could be ne- 
glected. However, recently, the slip length of PS on octadecyltrichlorosi- 
lane (OTS) and decyltrichlorosilane (DTS) coated Si wafers was discov- 
ered to range up to the scale of a micron 8, 9] 1101 111] . The dewetting 
patterns, in particular the shape of the dewetting rims around the grow- 
ing holes in the film were analyzed using a thin-film equation valid in the 
regime of large slip lengths, the so-called strong-slip model, [8] [12]. The 
thickness of the films was on the order of a few 100 nm and the dewet- 
ting mechanism was therefore nucleation rather than spinodal. However, 
since the hydrodynamic boundary conditions influence the rim shape it 
is to be expected that the power spectrum of spinodally unstable films 
is affected as well. Some dependencies of the dominant wavelength and 
time scale on the magnitude of the slip length for the case of an attracting 
van der Waals potential are discussed in [13]. In this study, we system- 
atically compare thin spinodally dewetting films with zero to large large 
slip lengths. Our main motivation is the availability of experimental sys- 
tems (PS on OTS or DTS-covered Si wafers U [10] [TT] ) which exhibit 
extremely large slip lengths and against which we can test our theoretical 
analysis in order to not only infer qualitative but also quantitative results. 
In particular, we consider an effective interface potential calculated from 
Hamaker constants as given in [7J. 

In the following Sec. [2] we first compare the dispersion relations for 
lubrication models for zero to moderate slip lengths with the dispersion 
relation for the regime of large slip lengths. We establish that even though 
the energetics is the same, films dewetting for large slip lengths have a 
qualitatively different dispersion relation as compared to sticky films, and 
therefore their structure factor has a maximum at a different wavenum- 
ber. We investigate the relevance of the difference between the dispersion 
relations for the sample systems PS on Si and on OTS(DTS)-covered Si 
in Sec. [3J and conclude in Sec. 2] 

For clarity of presentation we restrict our analysis to one-dimensional 
interfaces. The generalization to real two-dimensional interfaces is straight- 



forward for the lubrication models for zero to moderate slip. The general- 
ization of the lubrication model for large slip is not completely obvious due 
to the appearance of additional cross-terms and we have included it here 
in an appendix [14] . We note that omission of these cross-terms would be 
discovered at the level of the dispersion relation, where the growth rate 
would not only depend on the modulus of the wave vector but also on its 
direction, in contradiction to the isotropy of the physical situation. For all 
slip regimes we obtain the same dispersion relation for two-dimensional 
interfaces as for the corresponding problem with one-dimensional inter- 
faces, except that the wavenumber is now replaced by the absolute value 
of the wave vector. 



2 Spinodal dewetting 

2.1 The no-, weak-, and intermediate slip limit 

If the slip length b is small compared to the lateral length scale L in the 
dewetting film (i.e., the spinodal wavelength, see below), or comparable 
to L, the dynamics of a thin non- volatile Newtonian liquid film between a 
vapour of negligible viscosity and density, and an impermeable substrate 
is given in the lubrication approximation (i.e., for e = H/L <C 1, with 
the mean film thickness H ) by a degenerate parabolic partial differential 
equation of fourth order for the film thickness h(x,t) [3] 

dth = -d x {M(h)d x [Yl{h)+adlh]} (1) 

with the surface tension of the liquid- vapour interface a, and the disjoining 
pressure (DJP) H(h) — —dh$(h) (the negative derivative with respect to 
the film thickness of the effective interface potential) [TJ [2]. If the slip 
length is on the order of H or smaller, the so-called weak-slip regime, 
the mobility factor is given by M{h) — (ft 3 /3 + 6/i 2 ) /rj, with the fluid 
viscosity rj. For b ~ L, we have M(h) — bh 2 /r], called the intermediate- 
slip regime. The well-known no-slip regime is reached by taking the limit 
b — > in the weak-slip regime, leading to M(h) — /i 3 /(3 77), see [12] for 
more details. 

A homogeneous flat film of thickness H is linearly unstable if 8 2 uj(H) = 
—dhTi(H) < 0: In the early regime of dewetting, we can linearize Eq. {TJ 
about the base state H. For small perturbations Sh(x, t) = H — h(x, t) we 
get 

d t Sh = -M(H) dl [d h U(H) + o dlSh] . (2) 

The ansatz 5h(x,t) — 5h(q,t) exp(iqx) corresponds to a Fourier trans- 
formation with respect to x and leads to solutions of the form Sh(q, t) = 
<fi(g, 0) exp [co(q) t] with the dispersion relation 

u{k) = M{H) q 2 [O h n{H) - a q 2 ] 

= i (q/Q) 2 [2 - (q/Q) 2 ] . (3) 

For unstable films, i.e., for dhJl(H) > 0, there is a band of unstable 
modes with wave number < q < q c , with q c — ^[dhIl(H)/a. The 



dispersion relation has a maximum at Q — q c /s/2, which also defines the 
typical lateral length scale L = 2-k/Q which is also called the spinodal 
wavelength. The typical time scale, i.e., the inverse growth rate of the 
fastest growing mode, is given by T — 1/ [a M(H) Q 4 ] . 

The dispersion relation u(q) has the form given in Eq. Q in the no-slip 
regime, in the weak-slip, as well as in the intermediate-slip regime. For 
b < L changing the slip length therefore only changes the time scale T but 
not the position Q of the maximum of u>(q). Q, on the other hand, is only 
determined by the ratio of dhJl(H) and a. If one knows a, measuring 
the position of the maximum of co(q) for a number of film thicknesses 
allows to determine d}J\(h) and therefore the effective interface potential 
&(h) 6,7. This has been accomplished experimentally by measuring the 
power spectrum of the surface roughness S(q,t) — \Sh(q,t)\ 2 , which, in 
turn, can be calculated from the initial spectrum S(q, 0) = \5k{q, 0)| 2 and 
the dispersion relation Lo(q) 

S(q,t) = S(q,0)exp[2v(q)t]. (4) 

If S(q, 0) is flat in the range of unstable modes, then S(q, i) has a maximum 
at the same position as uj(q), i.e., at q — Q. 

2.2 Strong-slip limit 

For the case of a slip length b much larger than L, the thin film evolution 
can be captured by a different thin film model [121113] , called the strong- 
slip model by Munch et al. |12j . It can be written as 

■nu=4brjd x {hd x u)+bhd x [11(h) + d*h] 

— bh p (dm + ud x u) (5a) 

d t h = — d x (hu) , (5b) 

with the fluid mass density p and the horizontal flow velocity u(x,t). We 
note that this model is associated with plug flow in the cross- sect ion. For 
the experimental systems considered here we are not interested in the last 
term in Eq. (I5ap since the inertial term proportional to p is negligible. 
However, for completeness, we calculate the dispersion relation including 
this term. The first term on the right hand side of Eq. (I5a|l is proportional 
to the divergence of the total longitudinal shear stress component parallel 
to the substrate. Note that the velocity u cannot be eliminated from 
Eq. (0) even if the inertial term is neglected. 

If we perturb Eq. © about a resting [u(x, 0) = 0] flat film of thickness 
H we get to first order in the perturbation the problem 



r]5u=4bHr]dldu + bHd x [d h Ti(H) 9i + dl&i] 




— bH pdtSu 


(6a) 


dtSh — — H d x Su. 


(6b) 



With the normal modes ansatz 5h(q,i) = Sh(q, 0) exp [co(q) t + iqx] and 
Su(q,t) — Su(q,0) exp [cu(q) t + iqx] we get, after taking the derivative 



of Eq. ()6b|l with respect to x and subsequently eliminating Sh(q, 0), a 
quadratic equation for the dispersion relation uj{q) 



V uj(q) = - 4b H V q 2 u)(q) + b H 2 q 2 [d h U(H) - a q 2 



bHp[w(q)} 



with the two solutions 



(7) 



^ 2 ^ = ~j- p {^ 2 + m 



i±. i + 



Hq 2 [d h n(H)-aql] 



Uq 2 + 



(8) 



While the first solution u>i(q), corresponding to the plus sign in Eq. (JSJ, 
is negative for all q and does therefore not contribute significantly to the 
roughness spectrum, the second solution Ul2(q) is zero for q 2 [dfJI(H) — 
aq 2 ] = 0, i.e., for q = and for q = q c = \/2Q. For sufficiently small p 
or sufficiently large b the term under the square root is negative for large 
enough q > q c . The corresponding modes oscillate in time. However, 
there real part is given by Keu>2(q) — —V [4q 2 + 1/(6 H)] /(2p) < and 
thus, these modes are heavily damped. 

If inertia is negligible we can ignore the last term proportional to p in 
Eq. (or take the limit p — > in Eq. 0) and we get 



u(q) 



bH 2 q 2 [d h n(H)-aq 2 

77 (1 + AbHq 2 ) 
(q/Q) 2 [2 - (q/Q) 2 ] 
T[l + B(q/Q) 2 ] ' 



(9) 



with the time scale T = T)/ (ab H 2 Q 4 ) , as in the intermediate-slip regime. 
As expected, taking the limit B — > (b — > 0) we recover the dispersion 
relation of the intermediate-slip model. As compared to Eq. ([3]), we have 
one additional dimensionless parameter B = 4 & H Q 2 determining the 
relevance of slip, which is of order unity in the strong-slip limit discussed 
here. Both the time scale T as well as B depend on the slip length. 
Introducing the b- independent time scale T' = BT = 4r)/(a H Q 2 ) we 
can study the dependence of the dispersion relation on the slip length 
more easily. In the strong-slip model, neither T nor T' are the inverse of 
the maximum of the dispersion relation. Since the numerator in Eq. ((9} 
is proportional to the dispersion relation of the weak- and no-slip models 
discussed in the previous section and since the denominator is positive, 
the strong slip model has the same band of unstable modes q < q c . Taking 
the derivative of Eq. @ we get the position g max of the maximum of ui(q) 
at 



Qmax — ' 



yi + 25-i 

B 



(10) 



Fig. \T\ shows the dispersion relation for the strong-slip regime for various 
values of B. For B — > we recover the shape of the dispersion relation 
for the weak and intermediate slip regime, but the time scale T diverges 
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Figure 1: Dispersion relation in the strong slip limit in Eq. © (full lines) vor 
various values of B (see upper axis tics). The location of the maximum of the 
dispersion relation shifts to smaller values q/Q and the height approaches 2/T' 
from below for increasing B (dashed line). 



in this limit. For increasing B the location g max of the maximum moves 
to smaller values of q and the height of the maximum cj(g max ) approaches 
2/T' from below. In the limit B -> oo we get uj(q) -> [2 - (q/Q) 2 ] /T'. 



3 Experimental relevance 

In the experiments discussed in [B] the shift of Q with the film thickness 
was used to determine the effective interface potential &{z), assuming the 
dispersion relation in Eq. ([3]), i.e., for the weak/intermediate slip regime. 
The surface tension coefficient was a = 31 niN/m and the DJP had the 
form n(z) = -A/ (6ttz s ), with the Hamaker constant 2.2 x 1CT 20 Nm. 
From this we get for the dimensionless slip length B — 0.23 nm 2 b/H 3 . 
Therefore, in order to have, e.g., B > 1 or B > 0.1 for the lowest film 
thickness H = 2 nm in the experiment, 6 has to be larger than 35 nm and 
3.5 nm, respectively. 

The position of the peak in the power spectrum as a function of H 
and b normalized to the position in the weak/intermediate slip model that 
one would observe for the material combination studied in [6j is shown in 
Fig. [5] Cleary, the shift in the peak is larger for smaller film thicknesses 
H and larger slip length b. In order to get a deviation of the peak position 



on the order of 5% for the smallest film thicknesses of H — 2 nm, the slip 
length has to be larger than 8 nm, i.e., much larger than expected for PS 
on Si. 

If, on the other hand, the strong-slip regime was to apply and one tried 
to determine c^II from the measured peak position q ma x [see Eq. (|10|l ] with 
the equation valid only in the weak and intermediate slip regime, i.e., 
Qmax = \fdhU*] (2 a) with an "apparent" D JP II* , one would produce a 
systematic error in the measurment. The ratio of the actual DJP <9hII 
and the "apparent" DJP <9hIT is obtained by squaring Eq. (I10|l 



d h U*(H) = V1 + 2B 1 d h U(H). (11) 

In order to get a first estimate on the error we assume that only non- 
retarded dispersion forces are relevant and get to first order in B 

d h n*(H)^-^—- - f b u7 . (12) 

2 7r H 4 4 7r 2 a H 7 

Therefore, a spurious subleading term oc 1/z 5 is generated in the "appar- 
ent" effective interface potential &*(z) — — J z H*(z') dz' . 

The system considered in [6|, i.e., PS on a Si wafer covered with a 
native oxide layer, is known not to exhibit significant slip. However, re- 
cently it has been demonstrated, that covering the same wafer with an 
OTS or DTS brush leads to very large slip lengths up to the order of 
microns [HI [HI UHl- With the material parameters of OTS, SiO, and Si 
together with the thickness of the OTS layer and of the SiO layer in [8] 
we calculate the effective interface potential for a PS film on an OTS cov- 
ered Si waver using Eq. (3) in Ref. [7|. With this, we can calculate the 
deviation of the position of the maximum of the dispersion relation g m ax 
from the position in the weak-slip limit Q as shown in Fig. [3] With a film 
thickness of 4 nm a slip length of b — 1 /im is enough to generate a 20% 
shift in the maximum of the dispersion relation. Such a large shift should 
be detectable in the experiments. 



4 Conclusions 

In this paper we demonstrated that the hydrodynamic boundary condi- 
tion at the substrate surface significantly changes the power spectrum of 
film thickness variations in spinodal dewetting for experimentally relevant 
systems. Analysing only the peak position without knowledge of the hy- 
drodynamic boundary conditions can lead to significant systematic errors 
in the data analysis. As pointed out in |15| . viscoelastic thin films show 
a similar behaviour: while the position of the maximum of the dispersion 
relation is identical to the position in the Newtonian weak- /intermediate- 
slip case, in the strong-slip case it shifts to smaller wave numbers for 
increasing slip length 16, I5lll7|. 

The power spectrum of capillary waves should be affected by hydro- 
dynamic slip as well. However, up to now, a stochastic version of the 
thin-film equation is available only for substrates without slip [181 119] 
and the phenomenological ansatz taken in [20] can be extended directly 



'max 




H[nm] 4 






100 



80 



60 



40 



20 



/ / / 

/ / / 
/ / / 



/ 



/ 



i 



i ; i 



! ! i 

i i 



/ 



/ 



/ 



/ / / 



max 

/ 20% 

/ 10% 

/ 5% 

2% 

1% 

' 0.1% 



/ 




4 6 

fl[nm] 



10 



Figure 2: Shift of the position of the maximum g max of the dispersion relation 
ui{q) as compared to the position expected for the weak and intermediate slip 
model and experimental parameters from [6] as a function of the film thickness 
H and the slip length b. The contour lines indicate the slip length b necessary 
in order to obtain a deviation of 0.1%, 1%, 5%, 10%, and 20%. 
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Figure 3: Shift of the position of the maximum q max of the dispersion relation 
u>(q) as compared to the position expected for the weak and intermediate slip 
model and experimental parameters for a PS film on OTS from [7] as a function 
of the film thickness H and the slip length b. The contour lines indicate the slip 
length b necessary in order to obtain a deviation of 0.1%, 1%, 5%, 10%, and 
20%. 



to the case of weak- and intermediate-slip only. In the case of vanish- 
ing slip the position of the maximum of the power spectrum approaches 
Q from above as time proceeds |21j . Since the no-slip, the weak-slip, 
and the intermediate-slip case differ only in the mobility factor M, the 
same behaviour can be expected for the weak and intermediate-slip case. 
The mechanism for this noise-induced coarsening is simple: thermal fluc- 
tuations generate short wavelength fluctuations rather rapidly before the 
instability sets in, amplifying modes with larger wavelength. In the strong 
slip case, the maximum of the dispersion relation shifts to very large wave- 
lengths for increasing b, which could emphasize the effect of noise-induced 
coarsening. However, a detailed analysis of a stochastic strong-slip thin- 
film equation is needed to reach a conclusion on this point. 
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Appendix 



The generalization of the strong-slip model to 3D has been derived in |14| . 
With the two lateral velocity components u and v in the x and y-direction, 
respectively, the model equations in non-dimensional form are 

Re— =— [d x (4:hd x u + 2hd y v) + d y (hd x v + hd y u)\ 

+ d x [Ah + 11(h)} -^ (13a) 

Re— =— [d y (ihdyV + 2hd x u) + d x (hd x v + hd y u)\ 

+ dy [Ah + 11(h)] --^ (13b) 

dth = — d x (hu) — d y (hv) . (13c) 

Here, we abbreviate the total/materials derivative d/dt = dt+ud x + vd y 
and the two-dimensional Laplace operator A = d x + dy and Re is the 
Reynolds number. The lateral length scale L, the vertical length scale 
H, and the time scale T have been introduced in the main text. For 
this model the slip length is large and of order b — (3/e 2 , where /3 is an 
O(l) constant. The scale for the (disjoining) pressure is P — rj/T. The 
capillary number is Ca = T]L/(aT) = e. 

The linear stability of a flat film is again a straightforward calculation, 
by taking the first order in the perturbation, the problem for u — Su, v = 5u 
and h = H + Sh and making the normal modes ansatz (Su, Su, 5h)(q, t) = 
(&i(q,0),Su(q,0),Si(q,0)) exp[u)(q)t + iq ■ f], with q — (q x ,q y ) and r = 
(x,y). 



10 



Abbreviating Su(q, 0) = Suo, Sv(q,0) = Svo, Sh(q,0) — Sho, we obtain 
the linear eigenvalue problem 

77 Suo = — b f) H [(>iq x +q y )Su +3q x q y Sv ] 

+ ibHq x [d H H(H)-aq 2 ] Sho (14a) 

r/5vo — — brjH [(4q y + q' x ) Sv + 3q x q y Suo] 

+ ibHq y [d H n(H)-aq 2 ] Sh (14b) 

w(q) 5h = — iH(q x 5u + q y Sv ), (14c) 

with q 2 — |q| 2 . Here, we have neglected the contributions from the inertial 
terms and we have switched to dimensional quantities in order to connect 
to the main body of the article. From (|14p we find the same dispersion 
relation as we obtained for the 2D case in ([9]), i.e., w(q) = <±>(q), by setting 
the determinant of the matrix corresponding to the linear system to zero 
and by solving for uj(q). Alternatively, one can determine q x Suo + q y Suo 
from Eqs. (|14a|l and ()14b|l which yields 

{q x Suo + q y Sv ) (1 + 4bH q 2 ) — 

— q 2 [8 H U(H)-aq 2 ] Sh . (15) 

Inserting this in Eq. (|14cf) we recover Eq. ©. 
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